The bigravity black hole and its thermodynamics 
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We argue that the Isham-Storey exact solution to bigravity does not describe black holes because 
the horizon is a singular surface. However, this is not a generic property of bigravity, but a property 
of a particular potential. More general potentials do accept regular black holes. For regular black 
holes, we compute the total energy and thermodynamical parameters. Phase transitions occur for 
certain critical temperatures. We also find a novel region on phase space describing up to 4 allowed 
states for a given temperature. 



I. THE BIGRAVITY ACTION 



Bigravity is a theory of gravity with two independent dynamical metrics (denoted by g^^, and fp,u)- 
I ' The action, first considered by Isham-Salam-Strathdee [Ij], is 

1 /-r 1 

I[g^.u,f^.u] = / [V9R{9) + <^^fR{f)-U{g,f)\ (1) 

^ ■ where U{g, f) is an interaction potential depending both on g and /. The dimensionless parameter a 

measures the relative strengths of both Newton's constants. 

Bigravity has received intermittent but consistent attention since it was first presented. For some recent 
work see [2t6| . The first question one may ask is what is the physical metric determining the geometry 
of spacetime. This is the same as asking to what metric do particles couple to. We refer the reader to 
the literature for discussions on many different interpretations. In this work, we shall work entirely on 
vacuum space filled only by the two metrics. The problem we shall be concerned -regularity of the black 
hole horizon- is independent on the interpretation for the bigravity theory. 
The potential considered in [l| was ^ 



^ ; UigJ) = v,^{g^,. - Migo^p - UpWV"^ ~ FT^) (2) 

H ■ 

CZ • where ui has dimensions of mass. 

Flat space g^^i, = f/^^, = rj^^ is a solution to the equations of motion. The linear theory around this 
solution 

describes IJ a massless graviton^ h^i, + fxp^jy plus a massive one h^^i, — p^i, with a Pauli-Fierz mass term. 
The mass is equal to 

= i^^lii±^ (4) 



'Electronic address: jmaxbanados®fis.puc.cl] 
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^ In 1] y/J is replaced by which is consistent with diff invariance. We have set v = Q for simplicity, and to 

make contact with [8| . 
^ Note that a massless particle is expected for all choices of U because the action is diff invariant. 



2 



and the effective Newton's constant is G(l + a). Stability /unitarity of the hnear theory then requires the 
two conditions 

— >0, a>-l. (5) 

(T 

The potential ([2|) breaks the original diffxdiff symmetry down to the diagonal subgroup. The linear 
theory is nevertheless well-behaved because the mass term has the Pauli-Fierz form. [See [9.] for a general 
non-go theorem forbidding cross interactions between N gravitons preserving difF^.] 

The applications of bigravity to massive gravity has been extensively discussed in the literature and 
we shall omit here. Our main goal is to discuss the properties of black holes. Black hole thermodynamics 
in massive gravity (with a different action) has been studied in (To| . 

The main goal of this paper is to study in detail the prop erties of black holes solutions to the action 



([T]). We shall argue in Sec. |TT]that the solutions found in [ll|, associated to the potential ([2]), do not have 
regular horizons and thus they cannot represent black holes. There exists, however, a much larger class 
of potentials accepting exact solutions and having a unitary /stable linear theory [ij, In Sec. IIIII we 
consider a more general potential (with two parameters) and show that regular solutions do exist. We 
then compute the total energy, entropy, temperature and discuss the thermodynamical properties. We 
shall see that phase transitions exist when a < 0, which is allowed by unitarity/stability. 

II. THE ISHAM-STOREY EXACT SOLUTION 

Soon after the action ([T|) was proposed, solutions with spherical symmetry of the form 

g^^cix'^dx" = -h{r)df + ^ + r^dn^ (6) 

h{r) 

ff,^dx^'dx'' = ~X{r)dt^ + Y{r)dr'^ + 2H{r)dtdr + klr^dQ^ (7) 

were discovered by Isham and Storey'^ For our purposes here is enough to mention that h{r) and 

X{r) have the Schwarzschild-AdS form 

h(r) = 1 h Ar^ 

r 

X{r) = — (^1 - — + Ar2 j . (8) 

M, m and A are independent integration constants and fcp = |. The cosmological constants A and A are 
combinations of vi and A. See [ll| for full details. 

A remarkable property of the solution, already noticed in [llj, is that /^i, can also be brought to a 
Schwarzschild form via a coordinate change dt = dt' + n{r)dr, with a suitable ^{r). In the new coordinate 
system the metric f^^i, takes the form, 

A dr"^ 2 

U^dx'^dx'' = -X{r)dt'^ + —— + -r^dn^ (9) 

X{r) 3 

where A is a constant. The zero's of X{r) then represent horizons for the metric /^i/. The full solution 
is then simply two Schwarzschild metrics written in different coordinate system. 

The Isham-Storey solution fulfills another remarkable property. Consider an arbitrary linear combina- 
tion of both metrics 

qtLv = aQt^p + h f^,^. (10) 



^ Note that fcO 7^ 1 does not introduce conical singularities because these solutions are asymptotically (A)dS, and have 
horizons. We thank R. Mann for a useful conversation on this point. 
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Since /^jy is off-diagonal, the metric q^i/ is also off-diagonal. One can introduce a new time coordinate t' 
defined by dt = dt' + 57'(r)(ir, and again, for a suitable choice of ^'{r) the metric q^^ becomes diagonal. 
The remarkable property is that for any choice of a, b this metric again has the Schwarzschild form 

i q,^d.^d.^ - - (i - ^ + A'^') + i_2|r+Av + ''''''' ^''^ 

where A,B are constants that depend on A,B, and M',A' depend on all parameters M,m, A, a,b. 
This metric represents a black hole with a new mass parameters that depends on the particular linear 
combination. In particular, the location of the horizon depends on a, b. Thus, for fixed values of M, m, A, 
the Isham-Storey solution generates a whole family of black holes with horizons at arbitrary locations. 

Summarizing, at first sight, the Isham-Storey configuration has (at least) two horizons -one for each 
metric- defined by the points where the functions h(r) and X(r) vanish. These (candidate) horizons 
are located at different and independent points in spacetime because the zeroes of h{r) and X{r) are 
defined by independent integrations constants M and m. This interpretation is however not correct. The 
Isham-Storey configuration is singular because, as we now show, there is no coordinate system where 
both horizons can be made regular simultaneously. 

As a first attempt to prove regularity, the metric can be put in Eddington-Flinkestein coordinates 
making it regular at its own horizon. Remarkably, the metric /^i^ becomes regular at that point 
as well. However, this can only be achieved for either ingoing or outgoing coordinates (depending on 
the choice of sign for H{r)), but not for both simultaneously. As it is well-known only half of the 
Eddington-Flinkestein coordinates is not enough to declare regularity. 

Let us concentrate on the g^i/ horizon defined by the condition 

9M 

h{r,)^l +Ar2=0. (12) 

If A > this equation has more than one positive solution. The argument that follows applies to 
all regular horizons. We assume, however, that the horizon is non-extremal, that is h'{rg) ^ 0. The 

argument can be generalized to extreme horizons as well. Defining a proper radial coordinate dp — 
in a neighborhood of the horizon, the metric g^jy can be brought to the form 

ds^ ~ —a^p^dt^ + dp^ + angular part (13) 

where a is a constant. The 1-form dt is singular at the horizon (now located at p — > 0) but pdt is regular. 
The metric (jlSp represents a regular spacetime written in singular hyperbolic polar coordinates. 

Let us now look at the metric f^^, in the neighborhood of r^. The functions X, F, H take finite values 
at that point and thus, naively, the metric looks regular, 

d/^ = -X dt^ +ydr'^+ 2h dtdr + angular part. (14) 

Here x = X{rg), y — Y{rg) and h — H{rg) are non-zero for generic values of M,m and A. Then, 
since the 1-form dt is singular at r = (or /? = 0), the metric (jl4l) is actually singular there. A more 
geometrical way to see this is to note that the set of points (r = r^, — oo < t < oo) with fixed angular 
variables correspond to a timelike one-dimensional curve in the metric /^y, but a 0-dimensional point in 
the metric g^i.. This means that metric g^i, requires that the above set of points should be identified, 
while metric f^^i, require them to be all independent. The only way to make this metric regular is to 
assume that X{r) and H{r) also vanish at r = (with appropriate weights). A more explicit way to 
reach to the same conclusion is to put the metric (jl3p in regular Cartesian coordinates. This can always 
be done in a neighborhood of its horizon. Then one looks at the metric /^^ in this (regular) coordinate 
system. It follows that /^^ can be regular at the horizon of g^^i, if and only if its own horizon is also 
located there. 

In summary, if and r/ denote, respectively, the solutions to h{rg) — and X{rf) = we have 
concluded that 17^1, is singular at , and /^j^ is singular at Vg . The only truly regular solution is one with 
Tg = rf. This solution does exist, but corresponds to a proportional case /^^ ~ B g^^^ [B is a constant). 
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where g^i, is the Schwarzschild metric. This particular solution is the only regular black hole within the 
family of Isham- Storey solutions. 

It is instructive to recall the Reissner-Nordstrom black hole system having a non-zero Coulomb field 
Afj^dx^ — AqM. The same arguments exhibited above imply that the potential Aq^t) must be zero at the 
horizon if one demands the 1-form A^dx'^'' to be regular there. For the Reissner-Nordstrom black hole, 
this can always be achieved via a suitable gauge transformation. 

The Isham-Storey exact solution associated to the potential cannot be interpreted as a black hole, 
but only as the exterior solution to some mass distribution in bigravity. 



III. A BIGRAVITY BLACK HOLE 



We were motivated by the question of whether or not the singularity at the horizon was a generic 
property of bigravity, or a property of the particular Isham-Salam-Strathdee potential ([2]). The potential 
^ is one member within the infinite dimensional family of potentials ^ given rise to a unitary and 
stable linear Pauli-Fierz theory. It is then natural to ask whether or not other potentials may give rise to 
regular black holes. To answer this question (in a simple way without trying to prove a general theorem) 
we shall consider the 2-parameter family of potentials, 

U{g, f) = - ^^ly){9afi - fafi) x 

'MF^r^ ~ rr") + ^2{g^"9''^ - g^'g'^^)] ■ (15) 

Note that for all 1^2, this potential belongs to the Pauli-Fierz class as defined in [2i], and that for 1^2—0 
we recover the potential See [13] for a much larger class of potentials admitting exact solutions. We 
shall see that, for generic values of 1^2, regular, non-trivial (i.e., non-proportional) black hole solutions 
do exist. Since the condition 1^2 7^ represents an open set one can conjecture that, generically, the 
Pauli-Fierz family of potentials does accept regular black holes solutions. 

Before displaying the exact solution to the potential ([15]) we mention that (i) the potential ([T5l) also 
accepts the background 5^11/ = fp,i, — 77^,^. The linear spectrum contains a massless field plus a Pauli-Fierz 
massive theory. The mass is now. 



m 



(16) 



a 

As in the previous case, 1 + a must be positive to avoid ghosts. Then, linear unitarity /stability holds for 

> 0, a> 1. (17) 

a 

Our goal is to prove that a regular solutions exist, not to classify all solutions, which is a hard task. 
We consider the family of potentials where vi is positive while 1/2 is negative. We write them in the form, 

1/1 =p2, 1^2^ ~P'^. (18) 

The metric components are more complicated albeit exact [l3| . The metric ansatz is again given by 
©, and the functions f,X,Y,H are now given by the following expressions, 

h{r) = 1 Ar^ ^^—^ 

XW^ ' (l-5^-A,'=+2#3)^ ,19) 
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The exponent a and the cosmological constants are given by 



" = /o o,„2^/o^2 ^ (20) 



A 



2{3kl - 4fcg + 3) 
(3-2fc2)(3fc2-2) 
p'(fcg - l)(2fcgAt^ - 3fcgK^ + 3kl ~ 2) 



2crfc2(3- 2fc2) 



Here k is a new constant defined by 




P V - 2 



(21) 



Unlike the Isham-Storey solution, now the parameter fco entering in the ansatz ([T]) is an arbitrary constant. 

Note that this solution describes (anti)-de Sitter solutions even though there is no cosmological constant 
in the action. The main difference with respect to the Isham-Storey solution is the new term As 
shown in [isj . these terms are generic and appear for a large class of potentials. 

Given the functions h{r) and X{r) above, the function Y{r) is most easily expressed as the solution to 
the following simple relation 

X ^ 1 + klK^ 3(4 - 7fcg + 4fc4)afcg Q 

h k^K^ p^{2kl - 3)(3fc2 - 2)2 r"+2 ^ > 

Finally, knowing X(r) and Y{r), the function H{r) is given by 



^(^) = ±yfc2^-^(ono (23) 

These relations fix completely all unknown functions. We then have an exact solution to the equations 
of motion following from the action ([T|) with the potential (fT5|). 

Contrary to the Isham-Storey solution, the parameter k^ entering in the ansatz ([7]) is now an arbitrary 
integration constant. The range of this constant is severely restricted if we demand the charge Q to be 
localized. Since Q enters in the metric in the form ^ where a is given in ([20l) . Demanding a to be 
positive requires 

|<^o<y|. (24) 

A plot of q;(A;o) is displayed in Fig. ([T|). 

The cosmological constants appearing in our exact solution are not fundamental constants but depend 
on the state (see ([20]) ). We plot in Fig. [2] the values of A and A for the allowed range of fco given in (|24|) . 
We see that de Sitter and anti-de Sitter phase may coexist. We refer the reader to for a variational 
formulation on de Sitter spaces. Note that for tj < 0, which will turn out to be the most interesting case, 
both cosmological constants are negative for most values of fco, with a small de Sitter window. 

The main property of this solution is the appearance of 4 integration constants Af, to, Q and fco. One 
would like to give an interpretation to these constants in terms of asymptotic conserved charges. This is 
however not possible. As has been extensively discussed in the literature (see, for example, [ij), bigravity 
solutions break Lorentz (or de Sitter) invariance. In fact, in our case, the only asymptotic Killing vector 
is constant time translations. In this sense there is only one conserved charge at infinity. 

Interestingly, black holes do indeed have only one charge. As we see in next section regularity of the 
horizon imposes three restrictions on the values of Af, to, Q, feg- As a result, there is only one remaining 
free parameter that we can identify as the total energy (conjugate to time translations). 
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FIG. 2: Plot of A(fco) and A(A:o). Left panel a > 0; right panel a <0. 



IV. REGULARITY CONDITIONS 



The horizon for the metric g^v is given by the inner solution of the equation h{r) = 0. (The outer 
solution, if any, will be a cosmological horizon). As before, the metric f^i, has its own horizon defined 
by the equation X(r) — 0. As discussed before, the only way to have a regular horizon is to ensure that 
both X{r) and h{r) vanish at the same point. Moreover, the function H{r) -which appears in the metric 
multiplied by the singular 1-form dt- has to vanish at the horizon as well. 

There is one more restriction to achieve regularity. The "temperatures" of both metrics must be the 
same. This is a condition in either the Minkowskian or Euclidean sectors. Kruskal coordinates will 
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exists simultaneously for both metrics provided the derivatives /'(r) and X'{r) evaluated at r = are 
correlated. This is most easily seen in the Euclidean sector as a condition eliminating a conical singularity. 
But it is equally valid in the Minkowskian sector as a condition for the Kruskal extension to exist. 
Summarizing, let Vg the zero of h(r), 



Hrg) = 0. 



(25) 



Recall that h{r) is given in (|19|) . This equation merely defines rg without really imposing any restriction. 
Now, we demand the horizon of the f^,y metric to be located at the same spacetime point, that is 



Xirg) = 0. 

Next, as discussed above, we also demand the function H(r) to vanish at this point 

H{rg) = 0. 



(26) 



(27) 



Finally, the temperatures of both horizons must be the same so that a Kruskal extension exists for both 
simultaneously. This condition reads 

f{r)\r^^xoX'{r)\r^ (28) 

where xq is a complicated constant that we prefer to omit. 

Ideally, one would like to use (|25l26l27l32p to express, for example, M, m, Q, fco in terms of Vg. This is 
algebraically impossible. What we can do, it to use these equations to write M, m, Q as functions of fco- 
Thus, a regular black hole is built by choosing a value of ko, within the allowed range Once ko is 

given, all other charges 



M{ko),m{ko),Q{ko) 



(29) 



take definite values. The explicit formulas are analytical but not very illuminating so we skip them. It is 
more illustrative to present plots of A/(fco), m(fco), (9(fco) for the allowed range of ko. 

Fig. ([3]) shows M, m, Q for a positive and Fig. (HJ shows M, m, Q for a negative. These pictures are 
generic and do not change drastically for small variations of the parameters within the allowed ranges. 




FIG. 3: a = 0.5p= 1, p' = 0.5 
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FIG. 4: a = -0.7 p = 1, p' = 1.28 



One may be concerned that imposing all these conditions require severe fine tuning and the solution is 
highly non-generic. Note, however, that this fine tuning must not be imposed. The parameters Q, m, M 
of the solutions are not conserved quantities. Only one particular combination of them, the total energy, 
which will be discussed in the next section, is conserved during gravitational collapse. The particular way 
the collapse takes place will not be considered here, but it is a sensible guess, that in order to satisfy the 
cosmic censorship conjecture, the parameters will evolve to an equilibrium value that fulfill the constraint. 
From the canonical point of view we expect that this may also be derived from an extremum principle, 
in the same way that the constraint relating energy and temperature may be derived in the Euclidean 
treatment of black hole thermodynamics. We leave this to future work. 



We now turn to the problem of computing the mass, and thermodynamical functions for this black 
hole. 

Ideally, we should set up boundary conditions appropriated to the black hole problem, preserving some 
asymptotic group. However, bigravity solutions are neither asymptotically flat nor (anti-)de Sitter. For 
the black holes displayed above, each metric is asymptotically (anti-) de Sitter but with difi'erent speeds 
of light. There is no obvious asymptotic structure one would like to preserve. In order to move forward 
we do the calculation in the quickest way by allowing in the variational principle at least the family of 
metrics described by eh black hole, varying all parameters. 

Since the potential in the action has no derivatives of the fields, it does not contribute to the boundary 
terms. The total energy is clearly given by two ADM functionals, one for each metric. By usual methods 
p^l , we find for the variation of the Hamiltonian, 



V. TOTAL ENERGY AND BLACK HOLE TEMPERATURE 



5E = SE{g)+6E{f) 



= SM - 



5p'mk^{6k*-nk'^ + 6)dko 
p(3-2fc2)i/2(3fc2_ 2)3/2 



(30) 



plus terms that vanish on-shell. As explained above, the second line follows by plugging the black hole 
solution in the variation of the boundary term, varying all parameters. 6M is a expected contribution 
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from E{g) while the second term comes from E{f). 

It is very interesting to observe that both 5E{g) and 5E{f) have divergent contributions (as r — > oo). 
These are expected because the metrics have terms of the form Ar^ where A is not a fundamental constant. 
Instead, A depends on integration constants (see (1201) ) which are varied in the action principle. However, 
these divergent pieces nicely cancel each other and the sum SE{g) +dE{f) remains finite. A phenomenon 
like this one has been observed in other systems [l^. A general argument has been given in [l9|. 

Note that the second term depends on m and /cq. One may conclude that this solution has a new 
asymptotic symplectic pair ko,m (a charge with a corresponding chemical potential). This interpretation 
is however not correct. There is only one asymptotic time translation symmetry and only the total energy 
is a conserved charge. 

Integrability of the total energy then implies that m cannot be varied arbitrarily, but must be a function 
of kg. We have already seen that m is related to ko for regular black holes. We use this relation and 
write the total energy as 



V m(fco)fco2(6fc4-llfcg + 6 ) 

p A """^ (3-2fc2)i/2(3fc2_2)3/2 



where m(fco) is plotted in Figs. [3]and|4]as function of fco. We have arbitrarily chosen fco = 1 as a reference 
state. This is an allowed value and corresponds to the point where the exponent a{ko) takes its minimum 
value a{l) = 4. 

The temperature of the black hole can be computed either from the g^^, or /^^ black hole. Applying 
the usual formula for the (^^^ black hole, 



1 dh{r) 



47r dr 



(32) 



In Fig. [5] we plot the temperature for two choices of parameters a,p and p'. Focus first on the left 
panel with a = 0.5, p = l,p2 = 0.5. We see that for a given value of /3 (not too large) there exists two 
values of feg, and hence two values of M,m,Q. In the canonical ensemble, where /3 is fixed, there exists 
two black holes consistent with that temperature. To discern which one is realized one needs to compute 
the free energy. For negative values of a the situation is even more interesting. On the right panel we 
plot /3(fco) for a = —0.7,p = l,p2 = 1.405. For a given value of /3 there are up to 4 values of ko and thus 
4 states. Again, to discern which one is the most stable one we need to compute the free energy. The 
free energy will also tell us whether or not phase transitions among these states can occur or not. We do 
this calculation in the next section. 



VI. ENTROPY AND FREE ENERGY 

As shown in Fig. ([S]), for a given temperature there can be up to 4 values of /co- Hence, for a given 
temperature there may be up to four different black hole metrics. To decide which is thcrmodynamically 
favored we need to compute the free energy F — E — TS. 

The logic here mimics step by step the usual GR calculation, although the details are more involved. 
The full Hamiltonian action must be supplemented by two boundary terms at the horizon and infinity. 



/= / (bulk) + I3E-S{rg) (33) 

where E is given in (PT|) and (G = 1) 



S{rg) = TTrl+airklrl 

Aig) , A{f) 



(34) 
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FIG. 5: Left panel: a — 0.5, p = l,p2 — 0.5. In this case, there are two states (two values of fco) for each 
temperature. Right panel: a — —0.7,p — l,p2 = 1.405. Here there are up to four states for each temperature. 
This indicates possible phase transitions. The temperature T diverges near the extreme allowed values fco = 2/3 
and fco = 3/2- The horizon also vanishes at those points and the behavior resembles that of Schwarzschild. 



This formula for the horizon boundary term follows from general grounds (it can also be derived straight- 
forwardly). Indeed, the full action is a sum of Einstein-Hilbert actions. Therefore, the total entropy is 
expected to be a sum of the individual Bekenstein-Hawking entropies. Note that the second action is 
multiplied by u, hence the second entropy is multiplied by a. Up to this point, the horizon has no 
logical relation with the asymptotic charges M, m, etc. The action principle (|33|) is appropriated to black 
hole fields having a horizon at some arbitrary point r^. 

The next step is to evaluate the action on a saddle-point approximation. This has two implications. 
First, we solve the constraints plugging the regular black hole solution. The bulk pieces vanish (because 
there are a combination of constraints) and the total action becomes 

I^I3E-S{E) (35) 

where now the entropy does depend on the asymptotic charge. This is not the end yet. Since E is varied 
(/3 is fixed), there is one remaining equation SI/SE ~ which implies the first law: 



Of course, one can check that this value for /3 is consistent with the 'no-cone' condition 

The free energy — E — TS = F{I3) is a function of the temperature. Let us now evaluate F{f3) 
for various cases. As before, the formulas are analytic but not very illuminating. We display the main 
results with plots. Note that both F{ko) and /3(fco) (for all allowed kg) are known and hence we can plot 
F(l3) parametrically. 

We shall vary the couplings a,p and p' in a way consistent with stability/unitarity of the linear theory, 
as discussed in Sec. IIIIl For given values of a,p,p' we then vary ko within the range 

First, consider the case a = 0.5, p = l,p' = 0.5. Recall that for a positive p must be bigger than p' . 
The right panel in Fig. |6] shows F{/3). We recall that for a given value of (3 there are two values of ko (left 
panel) and thus two allowed states. As a consequence, the free energy is a multi-valued function of /3 and 
the lowest value determines the most stable state. In this case, the lowest F corresponds to the value of 
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kO 

FIG. 6: a = 0.5 and p ^ l,p' = 0.5 

/? with the biggest fco- This is true in the whole range of fcg and hence there are no phase transitions in 
this case. 

A more interesting situation occurs for a — —0.7,p — l,p' — 1.2 (for cr < 0, linear unitarity /stability 
requires p' > p). We plot in Fig. [7]F(/3). As before, there exists two states for a given (3. However, in 
this case, there exists a phase transition for (3 = (3* . For /3 > (3* , the most stable state is the one with 
biggest fco, while for (3 < f3* , the most stable state is the one with smallest fco. The physical free energy 
is the envelope of the curve F{(3) with a discontinuous derivative at /3 = /?*. Hence, a first order phase 
transition takes place at that temperature. 




FIG. 7: a = -0.7 and p=l,p' = 1.2 
An even richer situation arises if we increase the value oi p' , leaving p fixed. For a = —0.7,p = l,p' = 
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1.405 we plot in Fig. [8]the temperature and F{f3). For the range /3i < (3 < /32, four states are available for 
a given (3. The free energy F{p) (right panel) also exhibits this degeneracy and allows to determine the 
most stable state: for f3 < f3i and (i > P2 corresponds to the the biggest value of /cq, and for /3i < (3 < (32 
corresponds to the second biggest value of fco. 

In this case, there appears to be a "zero-order" phase transition when crossing the critical values j3 — /3i 
and (3 = (32- This hypothetical zero-order phase transitions has been discussed before in the context of 
self-gravitating gases (see [13] and [2l| for details). 




FIG. 8: cr = -0.7 and p = = 1.405 

Finally, in Fig. IHlwe plot the temperature and free energy for a sequence of values of p' , leaving p = 1 
fixed showing how the systems moves from two to four states. Since p' is a coupling constant in the 
action we would not interpret this eflFect as a second order transition. 



VII. CONCLUSIONS 



To conclude, we have argued that the Isham-Storey exact solution cannot represent a black hole but 
the exterior solution to a spherical star. This is not a generic property of bigravity, but of the potential 
chosen by Isham and Storey. For more generic potentials we find more integration constants in the 
solutions, which may be fine tuned to produce a regular horizon. We have computed the total energy, 
temperature and free energy of these solutions and have found several phases in the canonical formalism. 
Phase transitions do occur for certain critical values of the temperature, depending on the coupling a 
and the potential parameters p and p' . 
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